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. 5t , I consider quantum field theories that admit charged non-topological soli- 

tons of the Q-ball type, and use the fact that in a first-order cosmological 
phase transition, below the critical temperature, there is a value of the soliton 
charge above which the soliton becomes unstable and expands, converting 
space to the true vacuum, much like a critical bubble in the case of ordi- 
nary tunneling. Using a simple model for the production rate of Q-balls 
through charge accretion during a random walk out of equilibrium, I calcu- 
late the probability for the formation of critical charge solitons and estimate 
the amount of supercooling needed for the phase transition to be completed. 



1 Introduction 

A common feature of various cosmological models with a first-order phase 
transition is the generation of excessive entropy. This happens Because the 
probabihty of tunnehng from the false to the true vacuum is negligible when 
the supercooling parameter r] = {Tc—T)/Tc is very small. The question arises 
then if there exists some other mechanism that can facilitate the completion 
of the transition more rapidly, without the problems of excessive entropy 
production. 

In other cases of interest for cosmology a relatively large amount of su- 
percooling is actually needed in order to have out-of-equilibrium processes 
such as the generation of baryon asymmetry. A detailed picture of the phase 
transition is then necessary. 

Here I consider a different mechanism for a phase transition that is pos- 
sible in theories that admit non-topological solitons of the Q-ball type, and 
show that the rate of their nucleation can be large enough for the transition 
to proceed with much less supercooling than in other cases. 

There is a large amount of literature on non-topological solitons (see [|I| 
for a review). Here I will be interested in the case of Q-balls that are formed 
in the context of a first-order cosmological phase transition. 

The production of non-topological solitons in a cosmological phase tran- 
sition has been studied before P, P], and I will now describe briefly some 
previous results, focusing in the case of a first-order transition. 

Assume a potential U{(f), T) with a symmetry preserving minimum at 
= and a symmetry breaking minimum at </> = </>+, separated by a barrier 
with a local maximum Um at = (pM, and consider the evolution of the 
vacuum during the phase transition. At high enough temperature T the 
universe is in the symmetric vacuum. At temperatures close to the critical 
temperature T^, thermal fluctuations of the field become important. The 
fluctuating regions typically have volume V^ = ^^^, where ^ ~ Voir) ^^ ^^^ 
correlation length. The relative probability of finding a certain configuration 
of the field in a correlation volume is P ~ exp(— F/T), where F is the free 
energy of the given configuration. 

For high enough temperatures thermal fluctuations are strong enough to 
flip regions of the correlation volume from the symmetric to the asymmetric 
minimum and vice-versa. When the temperature falls below a certain tem- 
perature called the Ginsburg temperature Tg-, thermal fluctuations are not 
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strong enough to mediate such transitions. Since the number density of the 
particles is n^ ~ T^, we will have A^ ~ TqV^ charged particles in a correlation 
volume, and we expect a net charge fluctuation Qb ~ \/N in each bubble of 
the asymmetric phase (I assume a case with zero initial charge asymmetry). 

Usually, in theories that admit Q-balls, there is a constraint of a minimum 
charge Qmin that they have to exceed for reasons of stability. If Qb > Qmin 
then we expect a Q-ball to be formed in the bubble of the asymmetric phase. 

In the case of a first-order phase transition, when the symmetric vacuum 
becomes metastable, the potential energy density in the interior of the Q- 
balls becomes negative with respect to the potential energy density of the 
symmetric vacuum. Then if a Q-ball has large enough charge it becomes 
unstable and expands, converting space to the true vacuum, much like a 
critical bubble in the case of ordinary tunneling 0, |[ . 

In previous discussions of Q-ball production in the early universe, the 
constraint of the minimum charge, was quite restricting for the occurence of 
the phase transition, especially since most of the estimates were made using 
the thin-wall approximation for the solitons, which leads to a relatively large 
minimum charge [Q. 

It was recently shown that the thin- wall assumption is not correct, and 
that some theories may admit Q-balls without the restriction of a minimum 
charge, or with a minimum charge of order one. Then, charge accretion may 
lead to a growth in the number density of the critical charge solitons below 
some temperature, and thus mediate the phase transition f^. An estimation 
of the nucleation rate of these solitons is then important for the description 
of the phase transition. 

In the case where the solitons that are initially produced survive long 
enough to reach thermal equilibrium, and if a set of reactions that produce 
critical solitons remain in thermal equilibrium, then critical solitons are copi- 
ously produced. The conditions for this to happen are quite model-dependent 
[|, 1^ . Here I will deal with the case of critical charge Q-balls appearing before 
thermal equilibrium, and show that the rate of their production by charge 
accretion can be large enough for the phase transition to proceed. This case 
may be relevant also when we have soliton destruction or evaporation before 
they reach thermal equilibrium. 

I will begin by reviewing various approximations that already exist for 
the description of Q-balls in various ranges of the charge Q. Then I will 
focus on the case of a first-order phase transition in theories that admit 



non-topological solitons and calculate the rate of Q-ball production using 
results from subcritical bubble fluctuations and charge fluctuations. At tem- 
peratures near the critical temperature thermal fluctuations may flip regions 
from the symmetric to the asymmetric vacuum. If these regions also contain 
a charge fluctuation Q, and have energy less than Q free particles, then, 
instead of decaying into free particles, they will evolve into the soliton con- 
figurations. The evolution of these bubbles into solitons is affected by their 
interaction with the thermal environment. An absorption of a positive or 
negative charge can be simulated by a random walk with a + or — step. The 
difference from an unrestricted random walk is the existence of a minimum 
charge below which they disappear. This is the case of a random walk with 
an absorbing boundary. Using this model I will calculate the final soliton pro- 
duction rate. Once they are produced these solitons keep interacting with 
the medium, and their evolution can still be simulated by a random walk. 
Below the critical temperature there is also a critical soliton charge, above 
which the solitons become unstable. This is the case of a random walk with 
two absorbing boundaries. Using this model I will calculate the production 
rate of critical solitons. The critical solitons expand converting space to the 
true vacuum and mediating the phase transition. I will then estimate the 
amount of supercooling needed for the nucleation of a cosmological phase 
transition. The main result of this work is that, for generic values of the 
model's parameters, the phase transition can proceed much faster, with a 
much smaller amount of supercooling than other models (that do not admit 
nontopological solitons). 

This result may be important for cosmological phase transitions in theo- 
ries that admit non-topological solitons, such as the Minimal Supersymmetric 
Standard Model (MSSM) ||^, [0], or phenomenological theories of hadronic 
matter ||1|. 

In Sec. 2 I describe some approximations that already exist for the de- 
scription of Q-balls in various regimes [0, |^, 0. In Sec. 3 I will model the 
production of Q-balls during the phase transition, using results from subcrit- 
ical bubble fluctuations [§ and charge fluctuations 0. I will then calculate 
the probability for the formation of critical charge Q-balls through charge 
accretion and in Sec. 4 I will estimate the amount of supercooling needed for 
the completion of the phase transition. 



2 Non-topological solitons in various approx- 
imations 

Consider a quantum field theory of a complex scalar field with Lagrangian 

L = \{d,<pf-u{<p) (1) 

that is admits a global U{1) symmetry and a conserved charge 

0^e*V, Q = jJ(fx<p*dt<P (2) 

Minimizing the energy in the sector of field configurations that have a given 
charge Q yields the three-dimensional bounce equation 

v'0=f/:(0) (3) 

where 

f/.(0) = f/(0) - ^cuV' , <P{x, t) = e'^'(l){x) (4) 

This bounce equation has a solution when Ucj{(f)) has a nonzero global min- 
imum in addition to the local minimum at the origin. In order for this to 
happen, u must be less than M, where M^ = U"{0), and greater than some 
value c^o, for which f/(^(0) has two degenerate minima. Minimizing E^ with 
respect to u we get the soliton energy Eq. 

I will be interested in the case of nontopological solitons that appear 
near a first-order phase transition for which I will consider a generic effective 
potential depending on the temperature T: 

f/(0, T) = Im\T)<P' - A{T)<f + \<P' (5) 



where 
for which 



M^{T) = {aT^ - m^) (6) 

UM) = \{M^ - ^')0' - M' + \<P' . (7) 

The cubic term is understood as an effective term of the form (0*0)^' ^ that 
still has a global U{1) symmetry, and a is the square of the gauge coupling 



constant, assumed to be greater than A. The coefficient A will also depend 
on the temperature, but I will neglect this dependence near the critical tem- 
perature. 

In the limit of large Q, which is equivalent to u ^ ujq + 0, the form of the 
Q-ball solution simplifies considerably since Ui^{(f)) is such that the thin 
wall approximation is valid for the bounce solution. 

In the opposite limit of cj — >■ M we see that the barrier in the function 
Uuj{4>), which provides the bounce solution, becomes very shallow until it dis- 
appears. In this approximation, neglecting the A0^ term, and after rescaling 

i 

^ IT- ' ^ = e ..., . (8) 
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one gets a minimum 
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Eq = QM{1 - -) (9) 

QA^ 1 
- <o (10) 
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where S^ = 4.85 is the bounce action in the dimensionless potential ^ip^ — ip^. 
The bounce radius in C, is of order one so the soliton radius is 

1 1 

Rs , (11) 

Now we consider the case where the previous potential is at, or slightly 
below, the critical temperature where it has two degenerate minima at and 

f/(0) = A02(0-0+)2 (12) 

We consider large charge Q-balls in the thin wall approximation, but now 
the surface terms are important since the two minima are degenerate, and 
Uq = since the original potential has already a bounce solution. 
The Q-ball energy in the thin wall limit is 

E{Q) = ^ + 4vri?^^i - eV (13) 

where 

e = \U{<p+)\ S, = r^ d<p^2U{<P) (14) 



An approximate expression for the radius if we neglect the potential term is 

As the temperature drops further from the critical temperature, the neg- 
ative potential volume term becomes important, and there exists a value of 
e such that the Q-ball energy does not have a minimum, and is unstable 
against expansion. Equivalently, for a value of e, there exists a critical value 
of the charge Q, such that for Q > Qc the Q-ball is unstable and expands 
converting space to the true vacuum. These critical solitons have a charge 
Qc and a radius Re which are obtained from the simultaneous solution of the 
equations J| = g = to be § 
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3 Nucleation of critical Q-balls 

In this Section I will show that critical charge Q-balls can be produced 
through charge accretion from smaller charge Q-balls during a random walk 
out of equilibrium, and estimate the rate for this process. 

I will use the generic finite temperature potential of the form (^). At the 
critical temperature Tc we have 

V = -M^(j)^ - A(f)'' + X(j)^ = \<p\<p - <p+f 

where (^+ is the second degenerate minimum, and accordingly the parameters 
of the potential satisfy the relations 

M^ = 2A0^ , A = 2A0+ , M^ = J^ I2\ (18) 

The surface term is 



S, = J^ d<p pV{^) = ^ <P% (19) 



and at a temperature T = Tc — 6T slightly below the critical temperature 
the potential energy difference of the two minima is 

£ = \V{(f)+{T))\ = aT6T(f)l = aT^(j)lr] (20) 

where t] = ST/T is the supercooling parameter. 

I will estimate the production rate of Q-balls at the critical temperature 
using results of subcritical bubble calculations ^ and charge fluctuations 
[^]. Subcritical bubbles are bubbles of the symmetry breaking minimum that 
have smaller than the critical radius. In order to get an approximation for the 
production rate of subcritical bubbles we will use the ansatz for Gaussian- 
shaped bubbles of radius R 

Mr) = </>+e-^'/^' (21) 

and write an approximate expression for their thermal production rate 

r'{K) = aT^e-^^"^^)/^ (22) 

where a is an unknown numerical factor of order one and F{(f)f,) is the bubble 
free energy 

F(0,)= / 47rr^rfrU(^) + ^(0, T) | = 7/? + 5i?^ (23) 




For a bubble of the form ( pT]) and the potential (H) at the critical temperature, 
it turns out that the potential term will give the largest contribution to the 
production rate below, with a coefficient [§ 

6 = ^'1^1 {^M^ - y^A<p^ + ^0^) = c, <Pl M^ 

where ci = 0.26. 

Now, in a typical bubble of radius R there are 

N ^C2 T^R^ 

(j) and particles, where C2 = 3^2 = 0.51. The probability of a fluctuation 
of total charge Q in this number is [Q 



PiQ;N) 



N 



Then, changing variables from R to A^, the total production rate of bubbles 
of charge Q is 

r'(Q) = lTiN)PiQ;N)dN = 

= Syc^aT^^/^/g T ^?F^ /Rcii? (24) 

where from now on we write for the symmetry breaking minimum 0+. We 
can evaluate this integral using the saddle point method. The exponent has 
a minimum at a radius 

Rl = ^ . (25) 

Absorbing some numerical constants in the prefactor a of order one, and 
denoting a new numerical constant c = 2^/ci/^/c2 = 1.42 we get the final 
result 

The above expression for the nucleation rate T'{Q) of bubbles of the true 
vacuum with net charge Q is still not quite the same as the nucleation rate of 
Q-balls. First, the bubble energy should be smaller than the corresponding 
energy of Q free particles, otherwise the bubble will decay into free particles 
instead of evolving into a soliton. From our previous expressions we see that, 
for generic values of our parameters, this is indeed the case, and accordingly 
the charged bubbles, during their subsequent interaction with the particles 
of the heat bath, will preferably evolve to the soliton configuration instead of 
decaying into free particles. Although the initial bubble energy is not exactly 
the same as the corresponding soliton energy, the thermal fluctuations from 
the surrounding heat bath would rapidly compensate for that during the 
subsequent random walk. It should be noted, however, that the fact that the 
bubble energy is less than the corresponding energy of Q free particles is also 
a model-dependent condition and has to be checked in different models. 

Second, the Q-balls are bounces in a modified potential, and may be very 
different configurations from the true vacuum bubbles. If, for example, the 
bubble radius is much smaller than the corresponding soliton radius then it 
will expand until it becomes the lowest energy configuration of a Q-ball, but 
in the process it will meet a number of particles contained in this volume. 



and if there is a net charge fluctuation of negative charge Q, then the bubble 
charge will disappear. From the expressions for the corresponding radii of 
the subcritical bubbles and the solitons, (^) and (|TT|), and using (|TUp we get 
for the ratio of the volumes of the bubble and the soliton of charge Q: 

^'^' ^ e'M' = m'^ (27) 



Ms J ^c^cpMT A'^(t)T 

From (|l^) we have: 



thus 



with 



M^ ~ aT^ ~ A0% A ~ a0 (28) 



«'^'-40«',/? (29) 
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e ~ QA^/M"^ ~ Q\ (30) 



and for charges of order one we see that e is much less than one, in which 
case we indeed have the situation mentioned above. 

During the expansion of the bubble to a soliton, the bubble will encounter 
a number of positive and negative charges N which are contained in the 
volume of the soliton. This process can be simulated by a random walk, with 
a + (— ) step if a Q-ball absorbs a positive (negative) charge. The difference 
from an unrestricted random walk is the existence of a minimum charge 
below which the Q-balls are unstable and disappear. This is the situation 
in a random walk with an absorbing boundary, and the probability that an 



initial charge Q' will evolve to a charge Q is given by [10 



nQ',Q)='§,e-^ (31) 

where N is the number of particles contained in the soliton volume (from now 
on we will use Q and Q' to denote the units of charge above the minimum 
charge Qmin)- 

Hence, the final expression for the nucleation rate of solitons of charge Q 

is /.oo 

V{Q)= / V\Q')P[Q',Q)dQ' (32) 

Jo 



Using the previous equations (pG]), (plD, (|Tl|), we get the final result for the 
nucleation rate of solitons of charge Q 

rW) = rn.^T7^Qe-^-« (33) 

where 

and again we absorbed numerical factors of order one in the prefactor a. 

Once the Q-balls are formed and start interacting with the positive and 
negative charges in the heat bath, their subsequent evolution, before they 
reach thermal equilibrium or disappear, can be simulated again by a random 
walk, with a + (— ) step if a Q-ball absorbs a positive (negative) charge. 
The difference from an unrestricted random walk besides the existence of 
a minimum charge below which the Q-balls are unstable and disappear, is 
the existence of a critical charge above which the Q-balls are unstable and 
expand converting space into the true vacuum. This is the situation in a 
random walk with two absorbing boundaries. 

We are interested in the probability for a Q-ball, starting from an ini- 
tial charge Qq which is arbitrary but small enough so that the rate for its 
nucleation is not negligible, to reach the critical charge Qc before it reaches 
the minimum charge, and we also need an estimate of the total time that is 
needed for this process. 

Let X{Q,N) be the probability that a point at initial position Q will 
have crossed Qc by N steps before crossing 0. Let p and q be the relative 
probabilities for up and down motion in a single step, with p + q = 1. For 
sufficiently large number of steps X{Q, N) becomes independent of A^. In 
this case it is easy to calculate this probability X{Q). Assuming there is no 
asymmetry between positive and negative charges, and p = q = 1/2, we getQ 

m 

X{Q) = ^ (35) 

and the nucleation rate for critical charge solitons is 

T{Qc) = J dQX{Q)T{Q) (36) 

"'^Note that we have shifted the charges by Qmm so the exact expression for ( pq ) would 
be X{Q) — 7j-zo — ~ t)ut since Qc ^ Qmin the above expression ( p5[ ) would suffice. 
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For applications to a realistic situation we are interested in a different 
Icind of random walk, one where the time step is variable, depending on the 
value of the Q-ball charge at any step, since different charges have different 
cross sections. If n(0) and v are the density and velocity of the particles 
in the thermal bath and (j{Q) is the cross section for their absorption by a 
Q-ball of charge Q, then the time step varies as 

As a first approximation for cr(Q) we can take 

a{Q) ~ Ar^Rl (38) 

where Rq is the radius of the Q-ball. During a random walk from a small to 
a critical Q-ball this radius varies by several orders of magnitude, and also 



with different dependence on the charge, according to relations such as ([Tl|) , 

However, the simple result for the probability ( ^Sf ) is derived in the limit 
of a large number of steps N and turns out to be finite and independent of 
A^. It is the same whether we use the smallest or the largest time step, as 
long as the number of steps is large enough, and so it must be the same even 
if the time step is varying. 

I have performed several numerical trials of random walks with initial 
charge Qo and critical charge Qc, and time steps with various dependences 
on the charge and indeed it turns out that for any sample of random walks 
there is always a fraction of Qo/Qc of them that cross the critical charge, 
the dependence on time is very similar to the dependence on the number of 
steps, and they are all completed within Q^ steps. Some examples of these 
random walks are shown in Fig. ^ and Fig. 0. 

In these examples I used the numerical values of Qo = 10, Qc = 300, 
and for a sample of 3000 trials of random walks, a fraction of Qq/Qc where 
successful (crossed the critical charge). They where all completed within Ql 
steps and I show the results in Fig. |l|. I used a radius-charge dependence 
Rq oc Q^'^, for which the time step varies as At oc Q^^ , and I show the 
results for the duration of these random walks in Fig. ^ 

The assumption of a geometric cross section as an estimate for a{Q) is 



justified by recent results |^, at least for large values of the charge. From our 
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numerical trials of random walks we see that indeed the successful random 
walks spend most of the time in large values of the charge, so we expect it 
to be a good approximation. 

Another approximation that was implicit in the random walk model was 
that the absorption of a positive (negative) charge leads instantaneously to 
a Q-ball of charge Q ± I. Actually, it leads to an excited Q-ball of charge 
Q ± 1 according to reactions like 

(Q) + ^ (Q + 1) + AE{Q) . (39) 

However, one can easily see that generally AE{Q) is less than M, the el- 
ementary mass of the particles, and using (|37|), (^) one can see that the 
thermal fluctuations coming from collisions with the surrounding particles 
of the heat bath during an elementary time step of the random walk are of 
order T which is generally much greater than M, hence these collisions bring 
rapidly an excited Q-ball to its ground state. 



4 Nucleation of cosmological phase transitions 

From the previous expressions for the critical charge and the critical radius, 
Eqs. (p!6D, (0), and using the relations for our general finite temperature 
potential slightly below the critical temperature, Eqs. (|18|), (p!9|), (pOD, we 

with numerical constants C3 = 0.3, C4 = 0.7. For the parameters of the poten- 
tial M, A, 0+, I will use their values at the critical temperature, assuming as 
will be shown later that the transition happens at a small value of the super- 
cooling parameter. For numerical results I will take the coupling constants 
A = 0.1, a = 0.2, and the scale of the phase transition at the temperature 
T, = ITeV. 

I will assume a Friedman- Robertson- Walker (FRW), radiation-dominated 
universe, with the time-temperature relation 
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where 



/ 45 



and N = Nh + |A^/ are the essentially massless degrees of freedom. Thus, 
for temperature T = T^ — 6T slightly below the critical temperature the 
supercooling parameter is 

I will estimate the value of the supercooling parameter at which there are 
enough critical Q-balls produced, at least one per Hubble volume, in order 
to induce the phase transition |jT3[. 

The critical charge is infinite at and above the critical temperature. There 
are subcritical bubbles and solitons produced at these temperatures but the 
probability for their survival below the critical temperature is exponentially 
small, since this evolution is equivalent to a random walk with one boundary 
and for this process we expect an initial charge Q to disappear after N ^ Q'^ 
collisions. 

Let t{Qc) be the time needed for the completion of the random walk for 
a critical charge Qc- Then if 

r{Qc) ^ « Qc (43) 

the variation of the critical charge during the random walk is irrelevant, and 
the number density of the critical Q-balls produced at some time is 

n{Qc, t) = jdQ dt' g^r(Q, t') (44) 

where Q/Qc is the probability that a critical Q-ball with charge Qc will be 
nucleated starting from a Q-ball with charge Q. r{Q,t') is the nucleation 
rate of Q-balls at time t', and for a fast transition I will approximate it with 
the nucleation rate T{Q) at the critical temperature (]26|). 

From the previous results for the random walk we can also get an estimate 
for t{Qc) 

TiQc)^QlAt (45) 
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where, for relativistic particles with density n(0) in the heat bath, and cross 
section for absorption a{Q) ~ AtiRq, the time step is 

For an order of magnitude estimate of t{Qc) I consider an average charge 
Q = Qcf^, with the radius-charge dependence given by (p!5D, since it is near 
the phase transition but far from the critical charge: 



Also we have 



where, from (P|) 



p2 ^ 


1 Q4/5 


'^Q - 


Al/5 02 


dQc 


dQc dr] 


dt 


drj dt 


drj 


j2 



(47) 



(48) 



(49) 
dt 2Mp ^ ^ 

Using (H), (H), (H), (P), (pD, (H), the condition for the validity of our 
approximation ( ^31) becomes: 

dt 
^ ^'=T3 (Qe/2)4/5 2Mp 2 7/ ^ ^^ ^ 






1 / /l# \" a2 

A \^/aT ) TMp 

With Mp ~ 10^^ GeV and for a phase transition at a scale T ~ </> ~ 1 TeV^, 
this condition can be satisfied for a small value of the supercooling parameter 
77. If this condition holds, then 



n(g„t) = jdQdt'-^T{Q,t') 
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Here I took the lower limit of the time integral at the critical temperature, 
assuming that the condition (^) is quickly satisfied. 
From the previous expression ( P5| ) for T{Q) we have 

/•oo T21/2 „<t>M ^ 

I Q nQ) dQ = a p^^^ e ^ -TT Qm.n (52) 



SO we get 



2-17/2 _ <^A/ 
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"(«"• " = " WMpTI ^-'^ '^"■" (#) ^ ^'^P ""^ (83) 

where, = 0+ is the asymmetric minimum at the critical temperature, and 
M = M{Tc). Thus, the condition for a large enough density of Q-balls to be 
produced inside a Hubble volume in order to mediate the phase transition: 



gives 

with the numerical constant c = 1.42 and a a numerical factor of order one. 

If Qmin is small enough for this condition to be satisfied, then the transi- 
tion proceeds at the corresponding value of t]. Again, because of the factors 
of Mp, for a late time phase transition at the scale of ITeV, this condition is 
satisfied for a very small value of the supercooling parameter, if the minimum 
charge is small enough. 

The nucleation of a sufficient number of critical Q-balls depends sensi- 
tively on the minimum charge. We see that for large enough Qmin (greater 
than 22 in our example) the phase transition cannot proceed through Q-ball 
nucleation, at least not in the range of our approximations, something which 
is consistent with our expectations and with previous results 0. 

As Qmin becomes smaller, there is a sufficient number of critical Q-balls 
produced at small values of the supercooling parameter. 

For example, for Qmin = 20, we see that there are enough critical Q-balls 
produced per Hubble volume to nucleate the phase transition at values of 
the supercooling parameter 77 = ^ ~ 10~^. 
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At these values of rj, if we consider the case of ordinary tunnehng we 
find that the tunnehng rate is much smaUer than the expansion rate of the 
Universe, and the transition does not happen until much later times and 
larger values of the supercooling parameter. Indeed, in our range of values of 
the supercooling parameter t] ~ 10~^, the nucleation rate of ordinary critical 



bubbles at high temperature is calculated in the thin-wall approximation []T2 
r ~ T^e~^, where 

and we see that, for rj ~ 10"^, and for our values of parameters a and A, the 
result for the ordinary tunneling rate is indeed exponentially smaller than 
the Hubble expansion rate. 

For a better description of the phase transition one needs to calculate 
quantities like p{t), the fraction of the coordinate volume of the universe 
that remains in the symmetric phase = at time t after the beginning 
of the phase transition, which in our case may start right after the time tc 
when the temperature has dropped to the critical temperature Tc. The usual 
formula given in WM can be easily generalized to take into account the initial 



critical soliton radius, and the velocity of the soliton expansion: 

An (Rcih 



UMt) ^ - I *.r(a.«.)i^'(*0f (^ + I it. ^) (57) 



Here, R is the scale factor of the Robert son- Walker metric. Re is the critical 
soliton radius, and v is the velocity of the critical soliton expansion. The 
nucleation rate of critical solitons T{Qc,ti) and their initial radius Rc{ti) 
are calculated from (PO), (|^. Not much is known about the velocity of 



expansion of critical solitons. We do not expect it to be close to the speed 
of light, since we are very close to the critical temperature and the amount 
of supercooling is very small. However, if we assume that the expansion is 
similar to the expansion of a true vacuum bubble in the case of an ordinary 
first-order phase transition, then we can use a simple model to approximate 



it by balancing the pressure inside and outside the bubble [12 



C5— = c^Y^ri (58) 



where C5 = SO/tt^. 
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Substituting from the previous results, and keeping the first-order correc- 
tion coming from the initial soliton radius, we get 



+ X"V"{^)'^M (59) 



Here, the first term in the brackets would be the final result if we had ne- 
glected the initial soliton radius. The second term is the first-order correction 
that comes from the critical soliton size. 

This result confirms our previous estimates for the nucleation of the phase 
transition. The mechanism is similar to inhomogeneous nucleation or nucle- 
ation in the presence of impurities |]T5|. The corrections coming from the 
critical soliton size are negligible for generic values of the couplings in our 
model, but for exceptional values of the couplings, or maybe different models, 
they would have to be taken into account. 

5 Conclusion 



The final results (^), (|55|), (p9|) depend quite sensitively on the minimum 
soliton charge allowed for stability. Thus the relaxed minimum bound on 
the charge described in Section 2 is important. As explained before, the 
minimum charge that the solitons have to exceed for reasons of stability is a 
model dependent parameter, and our results for the occurence or not of the 
phase transition can also be viewed as conditions on this charge in order for 
the transition to proceed via Q-ball nucleation. 

I used a simplified model for the soliton production based on charge fluc- 
tuations, simulating the absorption of positive or negative charges by a ran- 
dom walk and ignoring other processes like soliton evaporation that may be 
included in more complete models. I assumed that there is no asymmetry 
between the Q-ball interaction with positive and negative charges, and also 
that the absorption of a positive (negative) charge by the Q-ball turns it 
instantaneously into a soliton of charge Q ± 1. I also assumed a geometri- 
cal cross-section for absorption of elementary particles by the Q-ball, which 
should be a valid approximation for macroscopic Q-balls of large charge. For 
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small charge this may break down but since the Q-balls spent most of the 
time of the random walk in large values of the charge this should be a valid 
approximation. 

Although these approximations are made to simplify the analysis, I do 
not expect the final results for the minimum soliton charge and the contribu- 
tion to the phase transition to change significantly if appropriate corrections 
are included, since these modifications may only affect the pre-exponential 
factors in (^). 

The most important correction is probably concerned with the model of 
soliton production via fluctuations of true vacuum bubbles. This is not as 
well-tested as other mechanisms of soliton production like, for example, the 



Kibble mechanism |TB|, but there are indications that the subcritical bubble 
formalism is correct [§. 
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Figure 1: In the vertical axis we plot the number of steps of the random 
walk, and in the horizontal axis the successful trials that cross the critical 
charge. 
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Figure 2: In the vertical axis we plot the duration time of the random walk 
in units of 0^/T^ and in the horizontal axis the successful trials. 
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